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The Landau-Pomeranchuk-Migdal effects on photon emission from the quark gluon plasma have 
been studied as a function of photon mass, at a fixed temperature of the plasma. The integral 
equations for the transverse vector function (f(px)) and the longitudinal function (gi(px)) consisting 
of multiple scattering effects are solved by the self consistent iterations method and also by the 
variational method for the variable set {po, Qo, Q'^}, considering the bremsstrahlung and the aws 
processes. We define four new dynamical scaling variables, x'^, x^, x\^, x%, for bremsstrahlung and 
aws processes and analyse the transverse and longitudinal components as a function of {po, Qo, Q^}- 
We generalize the concept of photon emission function and we define four new emission functions 
for massive photon emission represented by g^, Qt, Sl, Sl- These have been constructed using 
the exact numerical solutions of the integral equations. These four emission functions have been 
parameterized by suitable simple empirical fits. In terms of these empirical emission functions, the 
virtual photon emission from quark gluon plasma reduces to one dimensional integrals that involve 
folding over the empirical jr^'J^ functions with appropriate quark distribution functions and the 
kinematic factors. Using this empirical emission functions, we calculated the imaginary part of the 
photon polarization tensor as a function of photon mass and energy. 
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Quark gluon plasma (QGP) state is expected to be 
formed in the relativistic heavy ion collisions. In order 
to identify the plasma or a de-confined state, one needs 
to study the physical processes in quark matter, that 
can distinctly and conclusively identify this state of mat- 
ter, such as parton energy loss leading to jet suppression 
mechanism. In this context, electromagnetic processes 
such as photons and dilepton emission are also consid- 
ered as important signals. Photons and dileptons are 
emitted at various stages during plasma evolution, for an 
overview one may see [l|, |^ and the references therein. 
In depth study of photon emission processes in quark- 
gluon plasma were reported 0, Q including processes 
also from hot hadron gas Q ■ In the formalism of hard 
thermal loops (HTL) effective theory, the processes of 
bremsstrahlung ^7^ and a crossed process of off-shell an- 
nihilation called aws 0,0 contribute to photon emission 
at the two loop level. These two processes contribute at 
the leading order 0{aas) owing to the collinear singu- 
larity that is regularized by the effective thermal masses. 
Higher loop multiple scatterings having a ladder topol- 
ogy also contribute at the same order as the one and 
two loop processes 0, These higher loop rescat- 

terings, each giving finite decoherent correction to the 
two loop processes, need to be resummed. This resum- 
mation effectively implements the Landau-Pomarenchk- 
Migdal (LPM) effects arising due to rescatter- 



ing of quarks in the medium during the photon formation 
time. This results in an integral equation of a trans- 
verse vector function for the real photons. The photon 
emission rates are suppressed owing to the LPM effects 
fu\ llll | . The LPM modification of the photon spectrum 
is important at very low photon energies or at high pho- 
ton energies. For example, the bremsstrahlung radia- 
tion is strongly suppressed by almost 80% at very low 
qo/T values, whereas the photon emission from aws falls 
strongly for higher qo/T values HH. Thus the rescat- 
tering corrections modify the two loop contributions for 
bremsstrahlung and aws processes at opposite ends of 
the photon spectrum for real photons. 



I. EMISSION OF REAL PHOTONS 

The photon production rates from bremsstrahlung and 
the aws processes considering the LPM effects are esti- 
mated by using Eq^ llJ. 
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5Rf(p^) in Eq^ is the real part of a transverse vector 
function which consists of the LPM effects due to multi- 
ple scatterings. This can be taken as transverse momen- 
tum vector (p_L ) times a scalar function of transverse mo- 
mentum p± . The sign " denotes the dimensionless quan- 
tities in units of Debye mass mo as defined in • The 
function p± ■ SRf (p_l) is determined by the collision kernel 
{C{c[±)) in terms of the AMY integral equation (Eq|21l 
and the energy function SE (EqO. This was solved in 
[m using the variational approach. We simplified this 
variational approach and extended to finite baryon den- 
sity case 15). Moreover, we reported that the complex 
LPM effects can be very well reproduced by introducing 
thephoton emission function g(x) of dynamical variable 
X [ig. In terms of the single variable function g{x) , the 
photon emission rates are estimated as shown in Eq0| A 
figure for g(x) function for real photons will be presented 
later (Fig0. 



II. EMISSION OF VIRTIAL PHOTONS 

Processes that contribute to virtual photon emission 
in QGP at aag order ,17.] and the higher order correc- 
tions were reported. The processes qq gj* , and 
qg — > 97* contribute to photon polarization tensor at the 
order of aag and appear as the one loop processes in HTL 
method. Further considering LPM effects, photon emis- 
sion from QGP as a function of photon mass was also 
reported [T^]- For the case of virtual photons, these mul- 
tiple scatterings modify the imaginary part of self energy 
as a function of photon energy and momentum both, ft 
has been shown that the rescatterings modify the self en- 
ergy at low Q^, more importantly around the tree level 
threshold, — 4Af^ 19]. The dilepton emission rates 
are estimated in terms of the imaginary part of retarded 
photon polarization tensor, given by llQ|. 
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For the case of virtual photon emission having small 
virtulaity, the transverse vector function is determined 
by the EatTUland the energy transfer 5E{[tpj\_,pq, qq, Q'^) 
given in Ea llll ^19^]. For virtual photons, the coupling of 



quarks to longitudinal mode must be considered. This 
results in a scalar function of p_L which is determined 
by an integral equation and collision kernel in Eall2l 
For simplicity, henceforth I will refer to Eq^jas AGMZ 
equation. An analytical form for the collision kernel was 
given in [2(1 ]. owing to which the photon emission rate 
calculations became easy. The 6E is an energy denom- 
inator and can be interpreted as inverse formation time 
of the photon. For the case of massive photon emis- 
sion, this energy denominator is modified by replacing 
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In above equations, — po + qo and {,g are actually 
functions of pQ,qQ,Q^ represented as, f{p±,PQ,qo,Q'^), 
g(p±,Po,qo,Q^)- The Eg 1 101 and Eq|21 are identical ex- 
cept for SE energy factor. Further, Eg 1101 and the EgtT^ 
are also similar exc ept fo r the left side of AGMZ equa- 
tion is a constant v^ponj]- Aurenche et. at, solved these 
equations, based on a very elegant method of impact pa- 
rameter representation 19]. As mentioned earlier, we 
have solved the AMY equation by the variational method 
and a new method called iterations method. These meth- 
ods are formulated in terms of tilde variables. Therefore, 
it is required to transform the two above equations to 
tilde quantities. Following the definition of [ll| we have. 



Px = ^ ; f(Pi) = ^f(p±) ; S~E{p^) = ^5E{p^). 
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One can see that the collision kernel are related by 
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In the equations below, we intro- 



duce the calligraphic symbol C which is is a linear oper- 
ator and this should not be confused with the collision 
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kernel C{p±). The integral equations can be put in the 
following linear operator form. 
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The linear operator at a given momentum pj_ gives the 
sum of the transverse vector function differences between 
the current momentum and all other momenta arising 
from scattering weighted by collision kernel. The above 
linear operator suggests that the corresponding tilde op- 
erator is related to untilded operator exactly as the SE 
operator. It can be easily verified that for the transverse 
equation this means, C{p±) = (r/rn|j)C(p^). In the 
AGMZ equation when SE transforms to tilde quantity, 
its T/m% factor has to be absorbed by the transforma- 
tion of g function. Therefore, if we set the g to transform 

2 

as g{p±) — -ifi-g{p±), the equation for the longitudinal 
part could be written as. 
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This implies that the g is larger than g by a factor to|j. 
Therefore, when the solution for this EqEl is obtained 
and integrated over / j^^^g{pj_), the result will be 

larger than the true result by exactly r7i|, factor. This 
problem was not present for the transverse part because 
p^.f(p^) is independent of any to|, factor arising from 
the tilde transformation. Therefore, to take this anomaly 
into account correctly, we will introduce factor for 

the longitudinal contribution while estimating the imag- 
inary part of photon polarization tensor. Further, in the 
tilde variables, the above equation is not dimensionless. 
Therefore, we will divide the above Eq^l by m^i in or- 
der to get the following equation, where absorbing I /mo 
factor, g is now re-defined. 
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transforms as g(p±) = 
function. Therefore, the 



In the 

anomalous factor is now In the present work, we 

have solved the above Eott^ by iterations and varia- 
tional method. The variational method for longitudinal 
part has been re-derived We could get good agree- 
ment of distributions form iterations method and vari- 
ational results and these details will be discussed else- 
where. Following these two methods, we obtained the 



p^ distributions for the bremsstrahlung and aws cases 
for both transverse and longitudinal components. We 
considered five photon energy values, for each energy ten 
photon momenta and ten quark momenta, a total of 500 
distributions for transverse and longitudinal components 
of bremsstrahlung. Similarly the distributions were ob- 
tained for aws case. All these were compared in both 
iterations and variational method. 



III. EMPIRICAL ANALYSIS OF THE 
SOLUTIONS OF AMY, AGMZ INTEGRAL 
EQUATIONS 

In all the calculations that follow, we have used two 
flavors, three colors and as=0.30 and T=1.0GeV. In 
Eas l2UI2^ we define four dimensionless variables used in 
the following work. Here one may notice that the vari- 
able is the inverse of x variable used in ,16j for real 
photon production. are defined in the Eas l24l251 

which are obtained by integrating the p_L.f (p_l), g(p_L) 
distributions ( in Ref.^^ there was a 2 factor extra in 
the definition of / as given by 2p_L.f(pj^) ). are 
the quantities required for imaginary part of polarization 
tensor and the rate equations. 
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In the above equations superscript 6, a represent 
bremsstrahlung or aws processes depending on the po 
value used. The subscripts T, L represent contributions 
from transverse (f(pj^)) or longitudinal ((7(pj^))parts. 
Using the iterations data, we first obtained the integrated 
values of p_L distributions as given by of Eas l24l251 
We define the emission functions 5^"^ in Eg 1261 which 
are functions of variables x^"^. These x^"^ variables are 
given in Eas l27l3ni The g^'^ are obtained from corre- 
sponding values by multiplying with Cr^\ coefficient 
factors which are given in Egs l31l3^ Therefore, using 
the obtained from iterations and variational data, 
we constructed all the four g^"^ in Eq[2ni It should be 

noted that the quantities x^r^\ and Cr^\ listed in Egs l27^ 
are not definitions. They are the results of a hard 
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FIG. 1: The integrated values It{x) versus dynamical vari- 
able X2 defined in the text. The symbols represent the in- 
tegrated values of px distributions of about 500 cases of 
{po,qo,Q^} values. The different colored symbols represent 
different photon energies. 



FIG. 2: The dimensionless emission function gxix) versus 
dynamical variable x'^ defined in the text. The symbols rep- 
resent the integrated values of px distributions of about 500 
cases of {po, qo,Q^} values and are transformed by suitable 
coefficient functions. The different colored symbols represent 
different photon energies. The red curve is an empirical fit 
given in the text. 



search for dynamical variables hidden in the solutions of 
AMY and AGMZ equations. 
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Figure n shows the I^ix) as a function of X2 variable. 
All the bremsstrahlung transverse mode data consisting 
of about 500 points are shown in the figure. The figure 



exhibits some trends as a function of values, but 

still does not bring out the underlying dynamical scale 
hidden in the data. When plotted with respect to any 
other variables, the representation is much worse and no 
trend emerges from such plots. In Figure[21 we show the 
same data transformed to g\:{x\) as a function of vari- 
able which is a combination of X2 and xi variables. The 
hidden scale is revealed in this plot, where data for all 
different po, goi values merged into a single curve. Col- 
ored symbols represent different photon energies. This is 
the transverse mode bremsstrahlung emission function. 
Importantly, the crosses in the figure represent low 
data for transverse mode of the aws process, intention- 
ally plotted over here and will be explained in the next 
figure. We have obtained an empirical fits to this data as 
given by g^{x)(empirical) formula in Eg 1351 The need 
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FIG. 3: The dimensionless emission function grix) versus 
dynamical variable defined in the text. The required are 
the transformation coefficients and empirical fit (red curve) 
are given in the text. 



FIG. 4: The dimensionless emission function g\{x) versus 
dynamical variable defined in the text. The required cj, are 
the transformation coefficients and empirical fit (red curve) 
are given in the text. 



for the empirical fit will become clear later. 
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Similar exercise for transverse aws is shown in Figure 13 
However, we notice that the x variable is different from 
the previous figure. Further, the c?p that transform the 
integrated values into g functions are very complex. We 
have parameterized this data by an empirical function 
given in Eg 1361 One important thing to notice is for low 
enough < i,OGeV (possibly the data for < 4M^), 
the data deviates from general trends. These low 
data is shown by crosses in the figure. It was noticed 
that the same data when transformed as required in 
FigI21 and plotted versus x^, they are very close to the 



bremsstrahlung function. This is shown by crosses in the 
previous figure (Fig[5J). Therefore, the transverse mode 
emission functions for the bremsstrahlung and for the low 

aws are the same and are given by (/^(x^) (empirical). 
The Figure 0] shows the results for bremsstrahlung lon- 
gitudinal component. The corresponding coefficient 
functions and the empirical fits represented by (x) are 
given in Eas l33l37l Similarly, the Figure |5l shows the re- 
sults for aws longitudinal component. The correspond- 
ing c1, g1{x) are given in Eas l34l3(^ 

In Figure IHl we show the photon emission function 
in terms of dynamical variable x versus 5(2;) defined and 
shown in Ref.[l^ for real photons. We reproduce that 
data in the Figure. represented by black symbols. The 
inverse of this a;- variable of [l^ coincides with the present 
xi value and is related to the variables of the present 
study by a;~^(of Ref.p^Wxj^ Iq^^o • This is the motiva- 
tion for choosing the combination of xi and X2 variables 
for transverse modes in the present study; as —> 0, 
the dynamical variables of the transverse parts go over to 
real photon dynamical variable. Such a need is not there 
for longitudinal components. Further, the g{x) defined in 
|16| is two times larger than the definition in the present 
work for transverse modes (see factor 2 in 2pj_.f(p^) in 
EqQand not in EqsEE3- Therefore the rescaled ^g{x) 



versus x {= \ q2^q) for real photons is also plotted in 
the figure by blue symbols. For comparison, the results 
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FIG. 5: The dimensionless emission function g1{x) versus 
dynamical variable x1 defined in the text. The required are 
the transformation coefBcients and empirical fit (red curve) 
are given in the text. 

of our present work for versus are also shown in the 
figure by red curve. This red curve is the same red curve 
shown in FigEl given by g{x) = 10.0/ {5 + 3y/x + x). 
It should be noted that in this figure, the variable for 
X — axis is different for all the three curves. Further, 
we explained earlier, that for < A.OGeV^, the trans- 
verse component of the aws data lies very close to the 
bremsstrahlung curve (see crosses in the figure I^J. From 
these results it is obvious that for small quark momenta 
and low virtuality of the photons emitted, the transverse 
part of the virtual photon (7*) emission function and real 
photon (7) emission functions are the same as given by, 
(/^(x^)(7*) « g{x){'y) and this x for real photons is sim- 
ply X = x^ |q2=o- Further, the other longitudinal com- 
ponents for virtual photons are all negligible. Therefore, 
for high photon energies and small photon virtualities, 
the virtual photon emission is identical to real photons 
emission in the sense that their emission functions are 
the same. 



IV. EMPIRICAL ESTIMATION OF PHOTON 
RETARDED POLARIZATION TENSOR 

The imaginary part of retarded photon polarization 
tensor can be calculated using the integrated values 
as shown in Eg 1401 In the previous section, we used the 



FIG. 6: The dimensionless emission function 17(2;) versus dy- 
namical variable x taken from [T^. represented by black sym- 
bols. The half value of this g{x) in the inverse scale (that is 
same as present xi variable) is shown by blue symbols. This 
inverse variable coincides with the present when 0. 
The results for transverse mode of the virtual photon case ver- 
sus x^ variable is shown by the red curve (same as in Fig|2j . 

results from variational and iterations methods to ob- 
tain the /^''J^ values by integrating the distributions. We 

transformed these into g^"^ functions shown in Figs l2l5l 
and we fitted these by empirical functions Eas l35l3Hl The 
need for fitting arises, for example, for any {po,qojQ^} 
values, we can generate the integrated values l!^°j^{x) as 
given by Eg 1391 without solving integral equations. The 
imaginary part of retarded photon polarization tensor (in 
units of and represented by 311 i?) is calculated, as in 
Eg 1411 Note that in this section, no data from variational 
method or iterations is used. The empirical fit functions 
provide the required results circumventing the need to 
solve the transverse and longitudinal integral equations. 
The imaginary part of Hr has ^ factor for longitudinal 
part and this diverges as q^. This trend is seen in 

the Figs l7l8l as increases^. In this Eq^] one should 
note the multiplying factor (80asK/37r), in the present 



In this region, the formulae based on the kinematic conditions 
used will not be valid, and that the divergence serves as a good 
signal reminding us of the non-validity of the formulae (as sug- 
gested to me by F. Gelis). 
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FIG. 7: OTii plotted as a function of Q^/T^ for a 
photon energy of 50GeV. The transverse coinponents of 
bremsstrahlung, aws, the insignificant contribution from lon- 
gitudinal parts, and the 2 — > 2 processes contribution using 
the EqEH are all shown. The purpose of the figure is to test 
our empirical method compared to the results of Fig.3 of Au- 
renche et. al..\l^ shown by red circles. 



work it is equal to ^. Importantly, one should note the 
factor in the longitudinal part. This is because the 
g is already rescaled in Eq|Tni In Eq|?Tlthe superscript i 
denotes {6, a} bremsstrahlung and annihilation depend- 
ing on the value of the integration variable po. All terms 
in this equation are calculated and plotted in Figd to 
show the contribution to polarization tensor. 
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FIG. 8: Imaginary part of retarded photon polarization ten- 
sor for five photon energies represented by dashed curves in 
various colors obtained from the empirical emission functions. 
After adding the 2^2 contribution from Eg 1421 the results 
are shown by colored solid curves. 



The contribution of oneloop processes to imaginary part 
of the photon polarization tensor is given by equation 
below (formula taken from 'l{^). We will use this formula 
to evaluate the contribution of one loop processes and 
add to our empirical result to obtain the total result for 
retarded photon polarization tensor. 
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Figure [3 shows the 311^ plotted as a function of Q^/T^ 
for a photon energy of 50GcV. This figure has been com- 
pared with the results of Fig.3 of [13 represented by 
red symbols (approximately) for testing our empirical ap- 
proach. It can be seen that the empirical approach has 
predicted rather remarkably well the the results of 
over the Q^/T^ in the range of 1 — 100. Further, the 
polarization tensor increases rapidly with /T^ beyond 
100.0 The bremsstrahlung is insignificant obviously be- 
cause the qo is very high. The transverse component of 
aws only contributes, with a small contribution from lon- 
gitudinal component above > 20GeV'^. The 2^2 
processes contribution is calculated using the Eq021 The 
modification of the imaginary part of photon polariza- 
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tion tensor due to multiple re-scatterings in the medium 
has been well discussed in The LPM effects only 

marginally increase the 311^ at low . Here 2^2 
processes also contribute. Importantly, the rescattering 
corrections smooth out the discontinuity at the threshold 
— 4Af^. In order to understand how 311^ behaves 
with photon energy, we have calculated this quantity us- 
ing empirical approach for five photon energies. These 
results are shown in Figure |S| The dashed curves are the 
results from present calculations. When 2^2 contribu- 
tions are added to this, we get the solid curves. 

V. CONCLUSION 

The photon emission rates from the quark gluon plasma 
have been studied as a function of photon mass, consider- 
ing LPM suppression effects at a fixed temperature. Self- 
consistent iterations method and the variational method 
have been used to solve the AMY and AGMZ integral 
equations. We obtained the 5Rf(p^), 5ftg(p^) distribu- 
tions as a function of photon mass, photon energy and 
quark momentum. The integrated values of these distri- 
butions for bremsstrahlung and aws show a very good 
scaling in terms of different dynamical variables for these 
two processes. We defined four new dynamical vari- 



ables 2 for transverse and longitudinal com- 

ponents of bremsstrahlung and aws mechanism. In ad- 
dition, we defined four new emission functions namely 
g^{x^),g^{x^),g'l{x\)g1{x1), that describe the photon 
emission for these four components. We have obtained 
empirical fits to these emission functions that have been 
constructed using the exact numerical calculations. In 
terms of the new emission functions, estimation of the 
imaginary part photon polarization tensor reduces to one 
dimensional integrals. Using this empirical approach, we 
have calculated the imaginary part of retarded photon 
polarization tensor as a function of photon energy and 
mass. These self energy graphs inlcude the LPM effects 
due to rescattering in the plasma medium. For small 
photon values and high photon energies, the virtual 
photon emission function is a single function coinciding 
with the real photon emission function. 
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